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We show that a strongly interacting chain of Majorana zero modes exhibits a supersymmetric quantum crit¬ 
ical point corresponding to the c - tricritical Ising model, which separates a critical phase in the Ising 
universality class from a supersymmetric massive phase. We verify our predictions with numerical density- 
matrix-renormalization-group computations and determine the consequences for tunnelling experiments. 

PACS numbers: 71.10.Fd, 73.20.-r, 11.30.Pb, 74.55.+v 


Supersymmetry (SUSY) is a theoretical framework, which 
predicts a fermionic superpartner for every bosonic elemen¬ 
tary particle, explaining long-standing puzzles (see, e.g.. 
Ref. yj] and the references therein). The experimental ver¬ 
ification of SUSY, however, has remained elusive. Another 
important milestone of physics is Majorana’s prediction of 
fermionic particles that are their own antiparticle. Neutri¬ 
nos were Majorana’s original candidate for these so-called 
Majorana fermions, but their experimental status remains un¬ 
clear 12, [1]. 

While both Majorana fermions and SUSY are yet to be ob¬ 
served in high energy physics, Majorana fermions are pre¬ 
dicted to emerge as collective excitations in many-body sys¬ 
tems of electrons Of particular interest are unpaired 

Majorana fermions, localized to topological defects such as 
vortices or domain walls, that occur at zero energy. These 
Majorana zero modes (MZMs) provide a promising candidate 
platform for topological quantum computing ||2l- There has 
been significant recent experimental progress toward the de¬ 
tection of the condensed-matter incarnations of MZMs 113- 


1811 . Strongly interacting Majoranas may serve as building 


blocks for no vel p hases of matter, which remain relatively un¬ 
explored 119 - 2^ . 

On the other hand, there have been few works on the real¬ 
ization of SUSY in condensed matter physics |23-31|. The 
canonical example of emergent SUSY in statistical physics is 
the tricritical Ising (TCI) model in (1 -H1) dimensions 1 53 34 ]. 
This model is is the second simplest unitary minimal confor¬ 
mal field theory (CFT) in (1 + 1) dimensions. It has central 
charge c - ^ (Ising model with c - ^ being the simplest). It 
is also the only such CFT that exhibits SUSY. 

Similar to the Ising model, the TCI CFT has two re¬ 


alizations (i) spin models such as the Blume-Capel 

model 13514371] ■ in which, all local operators are bosonic, and 


(ii) fermionic models, in which both fermionic and bosonic 
local operators are present. Systems with Majorana fermions 
as local degrees of freedom provide promisin g ca ndidates for 
realizing the fermionic models of TCI CFT 1381] . The local 
operator content of the TCI CFT, which determines the exper¬ 
imentally accessible correlation functions, is directly related 
to the finite-size spectrum with periodic (antiperiodic) bound¬ 
ary conditions for the spin (fermionic) model. 
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FIG. 1. (a) The model with nearest-neighbor hopping and interac¬ 
tions between four nearest Majoranas. (b) The phase diagram of 
Hamiltonian O for g = 1 as a function of t. 


In the TCI CFT, there is only one relevant operator allowed 
by symmetry (independent of the realization), which must be 
fine-tuned to zero. This relevant operator can destabilize the 
TCI CFT, causing a phase transition either to a doubly degen¬ 
erate gapped phase or the Ising CFT phase (depending on the 
sign of the corresponding coupling constant). This operator 
preserves the supersymmetry of the Hamiltonian 1321] . SUSY 
is thus preserved in the gapped phase but spontaneously bro¬ 
ken in the gapless Ising phase OQI] . While nonsupersymmetric 
irrelevant operators could perturb the supersymmetric spec¬ 
trum deep inside the gapped phase, in the vicinity of the crit¬ 
ical point, the gapped phase is supersymmetric, with possible 
experimental signatures in tunneling experiments. 

Finding experimental realizations of the TCI CFT is of 
great interest. Important progress was made recently by con¬ 
structing models of coupled bosonic degrees of freedom and 
fermionic Majorana modes, with the intuition that these lo¬ 
cal bosonic degrees of freedom can serve as superpartners to 
the Majoranas 0. Purely fermionic field theories can also 
give rise to the TCI CFT |38, ^]. Here we show that indeed 
the simplest lattice model of interacting Majorana fermions, 
which may be realized in a superconducting vortex lattice 
and has clear experimental signatures in scanning tunneling 
microscopy (STM) experiments, gives rise to the TCI CFT 
and emergent SUSY. Similar to Ref. ItoI] . our model has the 
advantage that the TCI critical point can be reached by tun¬ 
ing only one parameter (thanks to translation invariance of 
the model). The two models, however, differ in an impor¬ 
tant aspect. The model of Ref. 0 can be thought of as 
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a model of interacting Majoranas upon integrating out the 
bosons. However, at the critical point, these boson-mediated 
interactions have a long-range character, whereas our interac¬ 
tions are strictly local. 

We begin by writing the Hamiltonian 



tty q Pi 

92-^00 ^ - C-OkS) - »Oo - • 

Pi dj ay+1 


yjyj+yyj+^yj+i. 


( 1 ) 


where yj - y\ (with {yi,yj} = 25,^) is the annihilation 
(and creation) operator for a MZM at position j in a one¬ 
dimensional lattice. This model can describe a vortex lattice in 
a narrow strip of a two dimensional topological superconduc¬ 
tor lEsIl . Throughout this paper, we focus on g > 0, for which 
supersymmetric phases emerge (setting g - 1 without loss of 
generality) lld^l . Positive g corresponds to attractive interac¬ 
tions between the underlying Dirac fermions, which may in¬ 
deed appear in the presence of superconductivity. The model 
exhibits a rich and co mp l ex p hase diagram for g < 0 , which is 
discussed elsewhere 143114411 . 

Our main result is the phase diagram shown in Fig.[T] For 
|f| > tc, the system is described by the Ising CFT with central 
charge c - while for 0 < |f| < tc, we have a gapped phase 
with broken symmetry. At \t\ - tc the system realizes the 
c = TCI model. As we will see, C is extremely small 
(relative to the interaction strength g). However, the regime of 
strong interactions is accessible in experiments due to a chiral 
symmetry (at chemical potential p = 0 in vortex realization 
of Majoranas), which forbids hopping processes ityjy/ 1123 , 


45h . By tuning 9 , we can then make t arbitrarily small without 
changing the interactions. 

It is convenient for the analysis of the problem to break 
the translation invariance of the system and write a more gen¬ 
eral Hamiltonian 


H ^iti ^ ajPj + it2 + .gi ^ aiPjOj+iPj+i 

j J J 

+ g 2 ^ /^y Q'y+ 1/5j +1 a j+ 2 , 


( 2 ) 


where aj = y 2 j and Pj = 727+1 ■ Each pair of Majoranas can be 
written in terms of one Dirac fermion Cj - (aj + iPj)l2. The 
product of two Majoranas is then related to the occupation 
number of a Dirac fermion through iajPj - liij - 1, where 
Hj — cjc/. Note that the pairing of MZMs into Dirac fermions 
is rather arbitrary and we could have introduced another set of 
Dirac fermions dj = (J3j + iaj+\)l2. 

In the limit of of f ^ cx) {gjt —» 0), we see from Eq. ([T]i that 
the system is described by a free massless Majorana theory 
corresponding to the critical phase of the transverse field Ising 
model; 




H, ^ iv dxiyRdjcyR - yidj^yi), 


(3) 


where a - 2 {yR + yi), P - 2 ( 7 * - yp), and v is a 
velocity related to the renormalized hopping t. No mass 


FIG. 2. The symmetry broken states at strong coupling. The thick 
lines indicate pairs of Majorana operators combined to form Dirac 
operators. The empty circles symbolize that the resulting Dirac levels 
are empty. For both of the states shown, there is another degenerate 
state, where are all Dirac levels are filled. 


term im J dxjRyi is present if we have translation symme¬ 
try. To incorporate the interactions into the effective the¬ 
ory, we Taylor expand the Majorana fields and obtain Hg = 
-256g f dxyi {Oxyi) 7 r ( 5 ^ 7 *) -!-■••, where the dots indicate 
terms of third and higher order in derivatives. From simple 
power counting, we find that the perturbation above is irrele¬ 
vant in the renormalization-group sense This implies that the 
Ising critical phase should extend at least to a finite g/t. 1411 ] 
Due to the large value of g/tc, the strong coupling limit of 
the Hamiltonian provides a good qualitative understanding of 
the gapped phase. It can be understood in terms of the occupa¬ 
tion numbers of the c and d Dirac fermions as shown in Fig.|2] 
First, we consider the case of f = 0. A dominant positive gi 
igi) gives ferromagnetic [all empty or all occupied] states for 
the occupation number of the c (d) fermions. The phase tran¬ 
sition at f = 0 and gi = g 2 (between phases with dominant gi 
and g 2 ) is expected to he first order so (i) all of these ferromag¬ 
netic states are present in the ground state manifold of Hamil¬ 
tonian ([TJ for f = 0 and (ii) there is a gap to excitations. As dis¬ 
cussed in the Supplemental Material (SM) 14711 . the presence 
of a first-order transition for g\ - g 2 and f = 0 follows from 
the connection of our model to spin chains and in particu¬ 

lar a generalized Ising spin chains with multispin interactions, 
which bears a strong similarity to the 8 -state Potts model & 


5311 . Despite being a small perturbation, the hopping terms 


(depending on their sign) lift the degeneracy between the all- 
occupied and all-empty states leading to a doubly degenerate 
(instead of 4-fold) gapped phase for 0 < |f| <s; g. 

The analytical arguments for a gapped phase at strong cou¬ 
pling are not rigorous. However, if the strong coupling phase 
is a doubly degenerate gapped phase, and assuming there is 
only one phase transition between the Ising phase and this 
gapped phase, then the TCI CFT is the most natural theory of 
this phase transition 11471] . The c = 1/2 Ising (c = 0 gapped) 
phase can be thought of as the line of a 2 nd-order (1 st-order) 
transition at t\ - t 2 (gi = ga) in a regime dominated by 
hopping (interaction). Similar to the Ising model with va¬ 
cancies 1-5411 ■ the critical point, at which the 2nd-order transi¬ 
tion changes to 1 st-order transition, naturally corresponds (at 
the mean-field level) to the <h® Fandau-Ginzburg theory at its 
tricritical point, whose universality class is given by the TCI 
CFT. 

This scenario needs numerical verification. As discussed 
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below, we found that the picture is indeed correct and the 
value of tc in the phase diagram of Fig. □ is tjg = 0.00405. 
This small value in turn implies a gapped phase whose short¬ 
est correlation length is thousands of lattice sites. Establish¬ 
ing the nature of the phases and determining the value of tc 
is therefore exceedingly challenging with most numerical di¬ 
agnostics such as extrapolation of gaps and order parameters 
as well as entanglement entropy. Despite this, we found that 
universal ratios in the finite-size spectrum provide a powerful 
numerical diagnostic for determining the phase diagram even 
though the system sizes we are able to reach are significantly 
smaller than the correlation length of the gapped phase in our 
model. 

Our evidence for the fermionic TCI CFT is the excellent 
agreement between the theoretical predictions for several uni¬ 
versal ratios at the critical point shown in the table below, and 
the numerically computed values of these ratios. 
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7 
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3 

8 

35 
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The subscripts A and P respectively indicate antiperiodic 
(APBC) and periodic (PBC) boundary conditions, Eo and Ei 
represent the energy of the ground state and the first excited 
state in a given fermion-parity sector (denoted by the super¬ 
scripts even and odd), and eo is the thermodynamic-limit en¬ 
ergy density in the ground state. The fourth universal ratio we 
use provides direct access to central charge c through the gen¬ 
eral finite-size dependence of the ground-state energy (which 
is in the even parity sector for APBC): = eoL 


Ittv _ 

L 12’ 


where L is the length of the system, if^ . 

We briefly outline the derivation of the above results based 
on the relationship between the operator content and the finite- 
size spectrum of the two CFTs. The results can be obtained 
from the formalism developed in Refs. 57] as discussed 
in the SM lld^ . We start with the Ising model, for which the 
predictions can be verified exactly in a free-fermion model. 
The Ising CFT has three primary fields I (identity), cr (spin), 
and e (energy) with conformal dimensions h - h - 0,-^, 5 , 
respectively [a held with conformal dimension {h, Ji) has scal¬ 
ing correlators - (x-vty^^(x+vty^^]. An 

integer (half-odd-integer) conformal spin h-Ji corresponds to 
a bosonic (fermionic) excitation. Now in a fermionic theory, 
we have APBC in the imagi nary time direction, which im¬ 
plies that modular invariance 15811 can be most easily satisfied 
if we also impose APBC in the spatial direction. The analog 
of the fermionic model in the Ising case is the free-Majorana 
model of Eq. Q, which has the conformal towers (I, I), (I, e), 
(e, I), and (e, e) with APBC (due to modular invariance 11581] 
as shown in the SM 14711 1, while the analog of the spin model 
has only diagonal conformal towers with bosonic excitations: 
(I, I), (cr, cr), and (e, e). 

In the fermionic model with APBC, the ground state and 
the first excited state of the even parity sector have operator 


FIG. 3. Top: The detection of the critical point through one of the 
universal ratios. Bottom: The values of four universal gap ratios at 
the tricritical point r/g = 0.00405 as a function of L (for systems with 
2L Majoranas). The numerical results (data points) show excellent 
agreement with the CFT predictions (hlack lines at 7/2,7/24,3/8, 
and 35/8). The DMRG truncation errors are insignificant (up to 1000 
states were kept in the computations). 


content (I, I) and (e, e) respectively (note that the conformal 
spin vanishes implying even fermion parity). Similarly, the 
ground state in the odd parity sector is doubly degenerate with 
operator content (I, e) and (e, I). For the Ising CFT, we then 
obtain= 2 £i;(^ + ^) andE°‘*^-E™ = 2 j^( 0 +l), 
which lead to the first universal ratio shown in the table above. 


The TCI model has 6 primary fields I, e, e', e”, cr, and 


cr', with scaling dimensions h - h — 0 , 


re- 


X 3 3 ^ 2_ 

10’ 5’ 2’ 80’ 16 

spectively. Similar to the Ising case, we have E^''™ - E^''™ = 
2|i^(E + X) However, in this case, (I, e) does not appear 
in the modularly invariant conformal towers of the fermionic 
realization of TCI CFT (notice that a conformal spin of ^ 
is neither an integer nor half integer). Here, the spin model 
also has 6 diagonal conformal towers, while, as shown in the 
SM the fermion model has 8 such towers with APBC, which 
include (e, e'), corresponding to the ground state of the odd 
sector with APBC 14711 . This has a conformal spin | = 5 . 

We then find E^**^ - E™ = ^(1 + ^), leading to the first 
universal ratio in the table above. The spectrum with periodic 
boundary conditions is a bit more involved but can be simi¬ 
larly derived using CFT methods (see SM ll47[l l. 


We numerically computed the four gap ratios above with 
the density-matrix-renormalization-group (DMRG) method. 
For the Ising and TCI CFT, the ratios above exhibit remark¬ 
able independence from the system size for large enough sys¬ 
tems. On the other hand, in the gapped phase (c = 0), at least 

^odd _^even 

one of the above ratios, namely grows with system 

^A.l ^A.O 

size for large enough systems (it has a linear dependence on 
system size for f = 0 as shown in Fig.O. This gap ratio was 
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used to detect the value of at the tricritical point. It plateaus 
at 7/2 for tc and approaches the Ising value of 1 /2 for larger t. 
Having found the value of tc, we then tested all four gap ratios 
for larger systems (see the bottom panel) and found excellent 
agreement with the theoretical predictions as seen in Fig. |3] 

We finally discuss the experimental signatures of the TCI 
CFT. STM provides a powerful tool to probe local density of 
states. Tunneling into MZMs can effectively probe the criti¬ 
cal exponent of the fermionic Green’s function. We start with 
the Ising CFT Q. The nonvanishing fermionic correlators 
are <rL,fi(b.r),yL,«(0,0)) = where b is a positive 

infinitesimal number. The equal-time fermionic Green’s func¬ 
tion then decays as 1 jx. A closely related quantity is the tun¬ 
neling current from an STM tip into a MZM, which goes as 
/[ oc y, where V is the bias voltage. 

In the TCI case, on the other hand, the leading fermionic 
operator corresponds to (e, e') with {h,h) - (3/5,1/10), 
which gives (x{t,x)x{0,0)) = 


to equal-time Green’s functions, which decay as |x| The 
tunneling current then goes as OTTII 


/tci oc sign(y)|y|'^^^ (4) 


In Fig. m we show the scaling behavior of the equal-time 
Green’s function of our model computed for TCI (f = tc) and 
Ising (t - 10^ ~ oo), where the predicted exponents are easily 
observed (.g = 1 in both cases). The only relevant operator 
(that induces a transition to the gapped phase from the TCI 
critical point) is (e', e') with dimension 3/5 + 3/5 = 6/5. We 
then expect a gap in the symmetry-broken phase that scales as 
(tc - f)^^^ near the critical point ll47n . 

As mentioned before, an important property of the rele- 
yant (e', e') operator is that it is supersymmetric 11381 l4lLl57 l. 
Therefore the SUSY of the critical point should extend into 
the gapped phase at least in the yicinity of the critical point. In 
the gapped phase, the power-law dependence of the tunneling 
current on V changes to exponential dependence, from which 
the gap to the leading fermionic excitation can be extracted. 
SUSY implies that the leading bosonic excitation has the same 
gap as the leading fermionic one. It should be possible to 
experimentally determine this bosonic gap from Cooper-pair 
tunneling yia a superconducting tip or other bosonic probes 
such as coupling to photons or phonons. 

Considering the effect of disorder on the rich physics of 
interacting Majoranas adds a new dimension to the prob¬ 
lem if^ : a recent manuscript, which a ppe ared shortly after 
the present paper, examines the effects 14311 . Our theory ap¬ 
plies to a translationally inyariant system. Experimentally, it 
is common to form Abrikisoy yortex lattices with translation 
inyariance due to energetic reasons. Spontaneous dimeriza¬ 
tion of the yortex lattice might occur which would indeed gap 
the system and destroy the TCI point. Howeyer, eyen in that 
case, if the dimerization is weak, some signatures of the criti¬ 
cal point suryiye in a crossoyer regime. 

In summary, we studied the phase diagram of the simplest 
model of strongly interacting Majorana zero modes in one 



FIG. 4. The scaling behavior of the fermionic Green’s function (ob¬ 
tained with DMRG keeping 800 states) for the Ising and the TCI 
CFTs observed in our microscopic model. 


dimension for attractive underlying interactions (which may 
be realized in the presence of superconductivity). Supported 
by extensive numerical calculations, we demonstrated that our 
model provides one of the few examples of emergent space- 
time SUSY in condensed matter physics and the first lattice 
model realization of TCI SUSY in a purely fermionic sys¬ 
tem with local interactions (a different type of SUSY has been 
predicted in a lattice model with 6-fermion interactions 1251] 1. 
The vortex lattice experimental realization of our model fos¬ 
ters the thus far elusive observation of SUSY, with clear signa¬ 
tures in the behavior of the tunneling current into a Majorana 
mode. 
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AR) and China Scholarship Council (XZ). 
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Supplemental Material for “Emergent Supersymmetry from Strongly Interacting Majorana 

Fermions” 


CONNECTION TO SPIN CHAINS 

A standard Jordan-Wigner transformation = Itij - 1 and 
crj = brings the Hamiltonian (2) of the main text 

to the form of a spin chain 

j j j J 

( 1 ) 

which for = 0 corresponds to the anisotropic next-nearest- 
neighbor Ising (ANNNI) model iU. For t\ - t 2 - 0, the 
Hamiltonian ([T]i is dual to a transverse field Ising model with 
4-spin interactions 

^ = ( 2 ) 

j J 

in terms of domain-wall variable = cF.cF.^j. 

The Hamiltonian ^ is self-dual al gi - g 2 Hi. As ar¬ 
gued in Refs, fli, the ground state of Hg is 8-fold degener¬ 
ate for .gi = 0, while the gi term causes quantum fluctuations 
between these degenerate states. This bears a striking sim¬ 
ilarity to the 8-state Potts model, which strongly suggests a 
first-order phase transition at gi - g 2 fli. 


SIMILARITIES WITH THE BLUME-CAPEL MODEL 

As mentioned in the main text, the canonical example of 
the TCI CFT appears in the Blume-Capel model, which is a 
quantum spin-1 chain with the Hamiltonian 

Hbc^-Yj + 1) - yS^ij) - . (3) 

J 

Similar to the one-dimensional spin-| transverse-fleld Ising 
model, which corresponds to the two-dimensional classical 
Ising model at finite temperature, this spin chain corresponds 
to a two-dimensional classical Ising model with vacancies 
{Sz = 0), with b serving as a chemical potential for these va¬ 
cancies. For b = -oo, the Blume-Capel model maps to the 
transverse-flled Ising model and exhibits a 2nd-order phase 
transition between a ferromagnet and a paramagnet with the 
critical point given by the Ising CFT. For y = 0, the model is 
classical and has a Ist-order transition between a ferromagnet 
and the S^ij) - 0 state. The line of the 2nd-order transitions 
terminates at the TCI point, where a line of Ist-order transi¬ 
tions begins. 

The above Hamiltonian has two dimensionless parame¬ 
ters. We can write a generalization of our model, which also 
has two dimensionless parameters, by staggering the hopping 
terms [Eq. (2) of the main text with gi = g 2 = .? > 0]. The 
phase diagrams are essentially the same as seen in Fig. 1. In 


1.2 

0.8 

0.4 


0 0.3 0.6 0.9 0 0.004 0.008 0.012 

b b/g 

FIG. 1; Left: the phase diagram of the Blume-Capel model (3) from 
Ref. 0]. Right: the phase diagram of our model with staggered hop¬ 
ping [Eq.(2) of the main text] with gi = g 2 = g > 0. The solid blue 
(dashed green) line represents a 2nd (1st) order transition with the 
TCI point denoted by a red star. 




our model broken translation fi - f 2 4 0 immediately gaps the 
system. 


FINITE SIZE SPECTRUM FROM CFT 


The finite size spectrum of the Ising and TCI model with 
both APBC and PBC can be derived using modular invari¬ 
ance 101 ■ Consider an imaginary time Feynman path integral 
for the partition function at inverse temperature /3 for free rel¬ 
ativistic Majorana fermions. As is well known, the boundary 
conditions must be antiperiodic in the imaginary time direc¬ 
tion. Thus, when we also impose APBC in the space direction 
the partition function 

Z = ^ e-2-r,^/L^ with E„ = 2^ (4) 

n 

is modular invariant. Each state corresponds to a left and 
right-moving factor and the energy is a sum of left and right 
moving parts: 


c 

yn ^ X„ + x„, with x„ = + Xi + n, (5) 

for nonnegative integers n, constants x,-, which characterize 
each conformal tower of states, and the central charge c of 
the CET. Each excited state in left and right moving confor¬ 
mal towers can occur independently. We define the character 
corresponding to a given chiral conformal tower as 

OO 

n=0 

where are nonnegative integers, which account for possible 
degeneracies of excited states in each conformal tower. The 
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partition function then has the form 

Z(J3IL) = ^ NijXi(J3IL)xj(f3IL). (7) 

ij 

Here the Njj are nonnegative integer encoding how many 
times each possible product of left and right conformal tow¬ 
ers occurs. Modular invariance imposes a restriction on the 
multiplicities Nij ||^ 

Nij-Yj^ii'SjfNrr, ( 8 ) 

i'j' 

where S is the modular S-matrix. 

For the Ising CFT, there are three conformal towers 


space. Za,p can be determined by a modular transformation 
from ZpA, where 

Zp,A=tr(-l)V^", (15) 

with F denoting the fermion parity. Notice that the opera¬ 

tors (/, e) and (e, I) have conformal spin and are therefore 
fermionic. To obtain Z^a, we simply insert a minus sign for 
all states with odd fermion parity. Thus 

Zp,A ^xl+xl- 2T2T3 ■ (16) 

Since Za,p is related to Zp,A by 

- 2 XiXjSii'SjrNl!^, (17) 

’■jJ'J' 


{cr,e,I) « (1,2,3) 


(9) 


we conclude that 


with the corresponding values of x, given by (1/16,1/2,0). 
The Ising modular S-matrix is 


S = 


0 -1/V2 1/V2 ' 

-1/V2 1/2 1/2 

1/V2 1/2 1/2 , 


( 10 ) 


There are two realizations of the Ising CFT with different 
local operator content, a spin model and a fermionic model. 
In the spin model, all three primary fields are present as local 
operators. A solution of Eq. ® for the multiplicities is given 
by Nij = 6 ij, which leads to the following operator content: 

(1,1), (cr,cr), {e,e). (11) 


As spins are bosonic, the partition function with periodic 
boundary conditions corresponds to the local operator con¬ 
tent. For the fermionic model, the field cr, which corresponds 
to the Ising order parameter, does not appear in the operator 
content since this is nonlocal in terms of fermions. Eq. (I8]l 
then has a unique solution as eliminating the first row and col¬ 
umn corresponding to cr leads to the eigenvalue equation 


/ 1/2 1/2 
I 1/2 1/2 


1/2 1/2 
1/2 1/2 


|1V> - |iv>. 


( 12 ) 


The matrix on the left-hand side has only one eigenvector 
{N\ - (1,1,1,1) with eigenvalue 1, which gives the operator 
content 


(1,1), (I,e), (e,I), (e,e). (13) 


= 2^^ (18) 

We can now extend this to the TCI model |@]. Eor the TCI 
CET, there are 6 chiral conformal towers 

(e,e',e",a,cr',I) « (1,2,3,4,5,6), (19) 

with (1/10,3/5,3/2,3/80,7/16,0) being the corresponding 
X,. The modular S matrix of the TCI model is given by 


-Si 

-Si 

S2 

V2ii 

- V2i2 

S2 

-Si 

-Si 

S2 

-V2si 

V2s2 

S2 

S2 

S2 

Si 

- V2s2 

- V2si 

Si 

V2ii 

- V2ii 

- V2i2 

0 

0 

V2i2 

- V2i2 

V2i2 

- V2ii 

0 

0 

V2ii 

S2 

S2 

Si 

V2s2 

V2si 

Si 


where sj = sm(j7T/5). Once again, we have a spin model and 
a fermionic model. Similar to the Ising case, the multiplici¬ 
ties of the spin model are given by Nij - 6ij due to modular 
invariance. We can then write the operator content of the spin 
model as 

(1,1), (e,e), (e',e'), (e'',e”) , (cr,cr), (cr',cr'). 

local in spin and fermion models, Z 2 even local only in spin model, Z 2 odd 

( 20 ) 

Similar to the Ising case, the operator content of the fermionic 
TCI model can be obtained form modular invariance and re¬ 
quiring that the nonlocal (in fermions) operators cr and cr' do 
not appear. Again, this leads to a unique solution to Eq. ® 
with the following local operator content for the fermionic 
model: 


Eor fermions, the partition function with antiperiodic 
boundary conditions corresponds to the local operator content 
and can be written as 

Zaa = xl'^xl'2X2Xi- ( 14 ) 

The spectrum with PBC for the fermionic model is given by 
Za,p, where the first subscript refers to time and the second to 


(1,1), (€,e), (e',e'), (e", e") , (e, e'), (e', e), (/, e"), (e",/). 

local in spin and fermion models, Z 2 even local only in fermion model, Z 2 even 

( 21 ) 

As in the Ising case, the local operator content of the 
fermionic model corresponds to the modular invariant spec¬ 
trum with antiperiodic boundary conditions |0 : 

Za,A - (xi +X2)^ + (X3 +X6)^- 


( 22 ) 
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Noting that operators (e', e) and (e", I) respectively have half¬ 
odd integer conformal spins 3/5-1/10=1/2 and 3/2-0 = 
3/2 (and are thus fermionic), we can write = (xi -^ 2 )^ + 
0f3 “ upon the insertion of (-1)^. From Eq. (14), we 
can then write the partition function with periodic boundary 
conditions as 


Z^,P = 20^4+4^5)- ( 23 ) 

From the above equations, we can read off the complete spec¬ 
trum with APBC and PBC for both Ising and TCI models. 
Note that the states with odd/even fermion parity correspond 
to the fermonic/bosonic operators^,. 

For Ising model, the ground state with APBC has x, - Xj - 
0 and is in the even sector. The degenerate ground states with 
PBC have Xi + Xi - 1/16 -H 1/16 = 1/8 and thus have en¬ 
ergy higher by (2nv!L){ \!%) than the unique ground state with 
APBC. The first excited states with PBC are a degenerate pair 
with energy (27rv/F)(l/4). For TCI model, the unique ground 
state with APBC also has xi - x; = 0. With PBC, the degen¬ 
erate ground states have energy higher by (27rv/L)(3/40) than 
that with APBC, and the first excited states with PBC have 
energy (2;rv/L)(7/8). 


TUNNELING DENSITY OF STATES 


Modular invariance also implies a relationship between the 
scaling dimensions of operators determining the exponents on 
their power-law decaying correlation functions and the param¬ 
eters x,„ Xn in the finite size energies. To each state in the 
spectrum with APBC there is a corresponding operator Ox„,x„ 
with <6);t,„jE„(T,x)Ox,„x„ (0,0)) oc ■ We use this 

general relation to determine the scaling of tunneling current 
with voltage in scanning tunneling experiments. 

Now let’s consider the single fermion Green’s function. For 
the free Majorana model (Ising model), we have two chiral op¬ 
erators (e,/) and (/, e), corresponding to Majorana operators 
jft and ji, where - {a ± P)I2. The only correlations are 
then 


<rfi,z.(b3c),y«,z,(0,0)) 


i 

Anivt ± X + id) ’ 


(24) 


where b is a positive constant of order the lattice spacing (one). 
Note that at equal times 


(ajao) = (/^jj3o) = 0 

(ajjSo) = (jSjao) = (25) 

In] 

On the other hand, for / = 0 but unequal times, 

<cro(0ao(0)) = mmm) =(26) 

2;7r(vf + lo) 

In STM we are interested in the tunn eling current between 
the sample and the normal tip, given by OlOl] 

</) = - 2 e^THmGR{-eV) 


where the retarded Green’s function 

Gfi(w) = -/ bfc'"'<[ao(0<Ao(0,«o(0)-Aj(0)]) (28) 

Jo 

is given in terms of Majorana fields ao and ordinary fermions 
i/ro describing the tip. It is convenient to obtain this retarded 
Green’s function by analytic continuation from the Matsubara 
Green’s function Gr{u)) - GM{io)n 6 + ico) which conve¬ 
niently factorizes 

Gm{t) = -<7’ao(T)ao(0))(r^o(T)<Aj(0)), (29) 


where T indicates time ordering. The Green’s function for the 
normal region at zero temperature is Gn(t) « -Dm It with 
Dm the density of states. For the noninteracting limit of the 
Majorana model -{T q'o(t)Q'o( 0)) ~ -1 /27rvT thus obtaining 


Gm(t) ~ 


Dm 

2 nvT^ 


(30) 


Performing the Fourier transform followed by the analytical 
continuation we find \mGR{—eV) — {Dmlv)eV. This gives a 
tunneling current 

/ = -27rrV—eV. (31) 

V 


Now consider the TCI model. We expect that 


Xr -> C\_XR\ + C2XR2 (32) 


where XRi the two fermonic primary fields of dimension 
(x,x) - (3/5,1/10) and (3/2,0) respectively. The Ci are 
nonuniversal constants. The corresponding correlation func¬ 
tions are 


Crfii(6-«kRi(0,0)> 

(XRiit, x)xr2(0, 0)) 


/ 


2;r(vf — X + /b)[(vf + iS)^ — x^Y^^ 


i 

2n(vt — X + i6)^ 


(33) 


The Green’s functions for xlj are obtained by taking x 
—X. The xi term in x should dominate at long times and low 
frequencies, so at equal times 


{aj/3o} oc 


e(j) 


(34) 


Using the Matsubara Green’s function 


<7’ao(T)Q'o(0)> oc 


1 

T|Tp/5 ’ 


(35) 


and following the same set of steps as in the Ising case we find 


ImGfi(w) oc Csign(w)|w|^^^ sin(77r/10), (36) 


where C is a positive real constant which can be expressed 
in terms of the T functions. The tunneling current, therefore 
reads 


(27) 


I oc sign(ey)|ey|^^^. 


(37) 
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Another quantity that can be predicted is how the crossover 
scale, from tri-critical Ising behavior, depends on f - tc- The 
leading bosonic operator consistent with the spatial parity 
symmetry is {e',e') of dimension 6/5, Oe/s- The effective 
Hamiltonian density near tc 

^^'Hrc+Ait-QOeis (38) 

where A is a constant of dimension (length)®^^. Under renor¬ 
malization, this coupling constant scales as 
Thus we can define a crossover length scale 

^^\t-tc\-^'\ (39) 

For f < tc, in the gapped phase, we should expect the gap to 
scale as 

A (X (f, - tf/\ (40) 
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